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Abstract 



We construct a new family of infinite-dimensional quasi-graded Lie al- 
gebras on hyperelliptic curves. We show that constructed algebras possess 
infinite number of invariant functions and admit a decomposition into the di- 
l/"") , rect sum of two subalgebras. These two facts together enables one to use them 

to construct new integrable finite-dimensional hamiltonian systems. In such a 
way we find new integrable hamiltonian systems, which are direct higher rank 
generalizations of the integrable systems of Steklov-Liapunov, associated with 
the e(3) algebra and Steklov-Veselov associated with the so(4) algebra. 



1 Introduction 



> 

The main purpose of the present paper is to introduce new integrable systems of 
Euler-Arnold type on the finite-dimensional Lie algebras. It is known, that such 
systems admit Lax-pair representations: 

^ = [L(t),M(t)}, (1) 

where L and M are some matrix depending on the dynamical variables 

Almost in all cases Lax operator depend also on additional parameter — so 
called "spectral parameter" w: L = L(w). This dependence permits to construct 
large number of integrals of motion via the expansion of functions I h (w) = TrL(w) k 
in the series of the powers of the spectral parameter. Usually dependence of Lax 
operator on w is rational or elliptic. 

The natural approach to construction of new types of integrable systems is search 
for the solutions of Lax equations with other (more complicated) dependence on 
the spectral parameter. We will solve this problem for the case of L(w) with the 
hyperelliptic dependence on the spectral parameter. 
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Our approach is based on the usage of infinite-dimensional Lie algebras. It 
is known 0,0, || that group theoretical explanation of the integrability of Lax 
equations on finite-dimensional Lie algebras with the rational spectral parameters 
is based on the loop algebras, i.e. algebras of the form g ® P(A, A~ x ). In the works 
0>HL HI it was shown, that Lax equation with the elliptic spectral parameters 
on the algebra so(3) and some its extensions could be obtained from the infinite- 
dimensional Lie algebras of the special elliptic matrix-valued functions with the 
values in so(3) via Kostant-Adler scheme 0- [||]. 

We generalize construction of on the case of classical matrix algebras of higher 
rank and investigate corresponding finite-dimensional integrable systems. Growth 
of the rank of algebra requires automatic growth of the genus of the curve. As 
a result we obtain algebras of gl{n)- , so(n)- and sp(n)- valued functions on the 
hyperelliptic curves of genus g, where n — 2g + 2 or n — 2g + 1, or, to be more 
precise, on its double covering. Obtained algebras have many nice properties. They 
are quasi-graded. They possess central extensions. They possess infinite number of 
invariant functions. In the rational degeneration of the curve they coincides with the 
ordinary loop algebras. So they could be viewed as "hyperelliptic " generalizations 
of the loop algebras. 

There exists other higher genus generalization of the loop algebras — so called 
Krichiver-Novikov algebras of the matrix- valued holomorphic functions on the 
Riemanian surfaces with two punctured points. But Krichiver-Novikov algebras do 
not admit Kostant-Adler scheme jl2) . Hence this generalization can not be used for 



producing new integrable systems with the spectral parameter on the higher genus 
curves. 

Contrary to the Krichiver-Novikov algebras, our algebras admit Kostant-Adler 
scheme, and hence, could be used for constructing integrable systems. 

Using our algebra as hidden symmetry algebra we obtain new integrable systems 
on g © g that describes two interacting generalized rigid bodies. This system is 
direct higher rank generalization of integrable case of Steklov and Veselov |13| on 
so(4) = so(3) ©so(3). We also obtain new integrable systems on the semidirect sum 
g + g, that generalize integrable case of Steklov- Ljapunov on e(3) fL4| . 

Structure of the algebra provides a possibility to give two- dimensional gener- 
alization of constructed integrable systems via finite-gap extension method 0-0. 
We will return to this problem in our subsequent publications. 



2 Quasi-graded algebras on hyperelliptic curves 
2.1 Construction 

Hyperelliptic curve embedded in C". Let us consider in the space C" with 
the coordinates wi, w 2 , w n the following system of quadrics: 

w\ - w) = a j -a i) i, j = 1, n, (2) 
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where a; are arbitrary complex numbers. Rank of this system is n— 1, so substitution: 



wl = W — CLi 



solves these equations. Moreover if we put y = Y[ w i we obtain the equation of the 

i=i 

hyperelliptic curve H: 

n 

y* = l[( w -a t ). (3) 

i=i 

Hence equations ||] define embedding of the hyperelliptic curve 7i in the linear space 
C n . Variable w is a local parameter on the curve and Oj its branching points. 

Example. In the n = 3 case all these objects have well-known analytical description. 
Indeed in this case curve under the consideration is elliptic. Its uniformization is 
made by the Weierstrass p- function and its derivative: w — p(u), y = l/2p'{u). 

Functions u>, are expressed via Jacobi elliptic functions |TIJ: W\ = — — — -, w 2 

dn(u) cniu) 
w 3 



sn[u\ 



sn(u) sn{u) 

Classical Lie algebras. Let q denotes one of the classical matrix Lie algebras 
gl(n), so(n) and sp(n) over the field of the complex numbers. For the subsequent 
we will need special form of their bases. Let us explicitly construct them. Let 
Jjj G Mat(n, C) be a matrix defined as follows: 

Evidently, a basis in the algebra gl(n) could be built from the matrices Xij = Iij, 
i,j G 1, ...,n. The commutation relations in gl(n) will have the following standard 
form: 

The basis in the algebra so{n) could be chosen as: X^ = 1^ — i,j G 1, ...,n, 
with the following commutation relations: 

[Xi t j,Xk t i\ = 5kjX it i — 5ijXkj + 5jjXk t i — 5k t iXj t i, 

and "skew-symmetry" property X^ = —Xji. 

The basis in the algebra spin) consists of the matrices X^ = 1^ — tiejl^i-j, i,j G 
— n, —1, 1, ...,n, with the commutation relations: 

[Xij,X k j] = 5 k jX it i — 5ijX k j + €i€j(5j_iX k _i — 8 k _iX_j t i), 

and additional property = —€i6jX_j_i, where ej = signj. 
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Algebras on the curve. For the basic elements X^ of all three algebras gl(n), 
so(n) and spin) we introduce the following algebra-valued functions on the curve 
H, or to be more precise on its double covering: 

X£(w) = Xij ® WiWj, X~(w) = Xij ® w^WiWj. (4) 

Here we put W-i = Wi in the case of sp(n). We will need the following definition 

Definition. Infinite-dimensional Lie algebra g is called Z-quasi graded, if there exist 
such p,q G Z + that: 

i 

g = ^g j; such that [Si,Qj} C ^ g i+j+k . (5) 

j'GZ k=—p 

The following theorem holds true: 

Theorem 2.1 (i) Functions X^Aw) andX^Aw) generate infinite- dimensional Z quasi- 
graded Lie algebra q-h . 

(ii)There exists a decomposition: q-h = 0-^ + 0-^, where areg^ andg^ are generated 
by Xf-{w) and X^(w) correspondingly. 

Proof. Let us introduce the following algebra-valued functions on the double cov- 
ering of the curve 7i: 

X™ = X i:j <g> w n WiWj. 

To prove the theorem we will need the explicit form of their commutation relations: 



Proposition 2.1 Elements X^ form closed algebra with the following commutation 
relations: 



1) [Xip X%\ = S kj X^ m+1 - 5 u X^ m+1 + atSaXg™ - o,*^ for the gl(n) 



(6) 



M A ij) A i,'J = <% A ;z ~ d ii x kj + d Ji X ki ~OikX n 



ij J ™3 il u u yl kj ~ u Jt- yl ki ^ik^^j 

L fcj _ a j°kj^u + diOik^ji 



+ OiSuXtf™ - afaXp™ + ai5 ik X"+ m - a&tXg™ for the so(n) (7) 



■3 J [ A ij> A ij'J - d fcj A iZ - °U A kj + e i e j[°3-l A k-i -Oi-k^j-i 

+ (HSaXJg™ - ajSkjX^™ + a i e i e j (a i 6 i - k X»* n - a^-_,X£+ m ) for the sp(n). (8) 
Let us put 



Q H = Span c {X"\ m e Z}. 



4 



From the Proposition |2TT] follows that g^ is quasi-graded. Besides, let us define the 
following subspaces: 

g+ = Span c {X™, m G Z + U 0}, q h = Span c {X^ m , m G Z + }. 

From the Proposition [2.1| follows, that they are subalgebras. Taking into account 
that X£(w) = X?-, X~j(w) = X~^ and the Proposition pTfl , it is easy to see that 
elements Xf-{w) generate subalgebra g^ and elements X^(w) generate subalgebra 
g^. Hence X^(w) and X^(w) generate g-^. Theorem is proved. 

Example. Let g = so(3). In this case constructed algebra will coincide with the 
"even" subalgebra of the algebra of hidden symmetry of Landau- Lifschits equations 
0. Indeed, putting X k = e^Xj,-, we obtain the following commutation relations: 

T V n V m ~\ <r yn+m+l , yn+m 



Remark. From the Proposition 2.1 follows, that in the rational degeneration of 7^: 



2/ 2 = w n , 

i.e when a, = 0, we obtain, that Wi = y/w and, hence, g^ = g, where g is ordinary 
loop algebra. 

2.2 Central extension 

To have the full analogy with the loop algebra, and keeping in mind possible applica- 
tion to quantum integrable systems, in this subsection we define central extensions 



of the algebras g^. First we remind the following definition [] I 



Definition. Cocycle \ on the quasigraded algebra is called local if 

x(Xij, x h) = 0) for all \n + m\> K and for some K G Z+. 

Let (X|F) = c n trXY, where c n = 2n for gl(n), c n = (n — 2) for so{n), c n = 
(2n + 2) for sp(n), be a standard invariant form on the classical Lie algebras gl(n), 
so(n) or sp(n). Let us define on g^ following bilinear form: 

X (X(w),Y(w) = j(X{w)\ d ^^)dw, (9) 

7 

where 1-cycle 7 in the complex plane of variable w encircles point w = 0. 

Although the above bilinear form is standard in the theory of loop algebras, it is 
not straightforward fact that it determines the cocycle on g^ due to the fact, that 
Wi are not one- valued functions of w. We prove this in the following theorem: 
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Theorem 2.2 (i) Bilinear form X is skew-symmetric and satisfies the properties of 
cocycle, and hence, determines a central extension ofg-^: gVc — &n + by the 
following formula: 

[X(w) + ac, Y(w) + (3c\' = [X(w), Y(w)} + X {X(w), Y(w))c. 

(ii) Cocycle X is local and its values on the bases elements are calculated by the 
following formulas: 

X {X™,X n kl )= X mn {a l) a 3 ){X l3 \X kl ) ) where (10) 

X mn (ai, a 3 ) = (n + l)5 n+m+2 ,o - (n + l/2)(a; + a 3 )S n+m+lfi + naia 3 5 n+Tnfl . (11) 

Example. In the case of rational degeneration = we obtain the following cocycle 
on the loop algebra: 

X (X™, X£) = (n + l)S n+m+2 , (X i3 \X kl ). 
It pass to the standard one after renaming the indices: m + 1 — > m, n + 1 — > n. 

2.3 Coadjoint representation 

To define the coadjoint representation we have to define g^. We assume, that 
0^ C Q Cg> A, where A is an algebra of function on the double covering of the curve 
7i. Let us define invariant pairing between L{w) G q\ and X{w) e Qn in the 
following way: 

(X(w),L(w)) f = c n res w=0 f- 1 (w)y~ 1 (w)(X(w)\L(w)), (12) 

where f(w) is arbitrary function on the curve H. It is easy to show, that element 

dual to Xr. m with respect to this pairing is Y£ = (X^Y = Jy Jyy ' X*,. 

n *3 v 13 > w . w . u 

Hence the general element of the dual space has the following form: 

iw = E E i^ wm ' lf{w)y{w) x: 3 (is) 

v ' ^ ^ 13 w i w j 13 v ' 

m£Zi,3=l ■> 

Coadjoint action of algebra on its dual space g* u coincides with commutator: 

ad* x(w) L(w) = [L(w),X(w)}. (14) 
Explicitly coadjoint action of algebra has the following form: 

1) K, Y$\ = 5 k3 Yt m+1 - 5 u Y k n 3 +m+1 + a 3 5 u Y k n + m - a^Yt™ for the gl(n) 
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2) [X$,Y2] = 6 kJ Y2 



n+m+l 




n+m+l 




n+m+l 



- 6^ 



rn+m+l 



+ a 3 5 a Y£+ m - ai 5 kj Y^ m + a 3 5 lk Y; +m - a^Y^™ for the so(n) 

3) [X%, Y k 7) = 5 kj Yt m+l ~ S U Y k f m+1 + e,e^-_^_r +1 - ^ k Yf+ m+1 ) 

+ a 3 5 a Y k n+m - ai8 kj Y£ +m + e i e i (a^ i _ fc y/+ m - a^Y^™) for the sp{n). 



It is evident from the above formulas, that g^ is a quasi-graded g^-module: 



where (g^) m = S , pan {ly*|i, j = 1, ...,n}. 

Remark 1. Putting f(w) = y(w) in the definition of the elements of dual space we 
obtain, that g^ C g w . If /(it;) 7^ then generally speaking, spaces g^ and g-^, 

viewing as subspaces in the algebra of functions on the double covering of the curve 
TC, do not coincide. 

Remark 2. If / = 1 and g = so{n) elements of q-h + g^ form a closed algebra. This 
will be analogue of the algebra of hidden symmetry of Landau- Lifschits equations 
0. Unfortunately, for n > 4 it does not admit a Kostant-Adler scheme and can not 
be used for the construction of integrable systems. That is why we will not consider 
it here. 

From the explicit form of coadjoint action QI^D follows the next statement: 

Proposition 2.2 Functions l!^(L(w)) = res w= ow~ m ~ l TrL(w) h , where m eZ, are 
invariants of coadjoint representation. 

Hence constructed Lie algebras not only admit decomposition into the direct sum 
of two subalgebras but also possess infinite number of invariant functions. This 
permits to use them in construction of integrable systems. 

3 Integrable systems from hyperelliptic algebras 

3.1 Poisson structures and Poisson subspaces. 

First Lie-Poisson structure. In the space g^ it is possible to define many Lie- 
Poisson structures using different pairings. We will use the pairing ([12]) with f(w) = 
w: 




(X(w),L(w))-i = c n res w=0 w 1 y 1 (w)(X(w)\L(w)). 
It defines brackets on -P(g^) in the following way: 



(15) 



n 



dG OF 




-1 




(16) 



l,m£Z i,j,p,s=l 



The following Corollary of the Proposition |2]2| holds true: 
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Proposition 3.1 Functions 7^(L(u>)) are central for brackets { , }. 

Taking into account that = (L(w), X^ m )_i, it is easy to show, that for the 
coordinate functions /™ these brackets will have the following form: 

1) 05. *y } = « +m_1 - <M# m_1 + MiiC - «AA" +m for the ff Z(n) (17) 

2) (£} = k 3 lT m ~ l - Sul^™- 1 + SjJ^™- 1 - S^' 1 

+ a t 5 t il^ m - a 3 5 k] ll +m + a t 5 tk l]? m - a^l n k + m for the so(n) (18) 

o\ fin rrn\ _ r m+m-1 r 771+m-l , / r m+m-1 r m+m-l\ 

+ aMkt m - a Ajlu™ + ^A^-kl]T - for the sp(n). (19) 

From the explicit form of the Poisson brackets follows, that Lie-Poisson brackets 

oo oo 

in the subspaces (g^)_ = Yl (dn)-m and (g^)+ = (dh)m, up to the reverse 

m=0 m=l 

of the sign of the upper indices, repeat commutation relations of the algebras g^ 
and g^. correspondingly. This, of course, is the result of the following duality: 

(&)• = (3h)-. (0w)* = (&)+■ 

Second Lie-Poisson structure. Let us introduce into the space g^ a new Poisson 
brackets { , }o, which are a Lie-Poisson brackets for the algebra g^, where g^ = 
g w 0g^. Explicitly, this brackets have the following form: 

Oo = -{iS-.O, n,m > 0, {Z?.,Oo = O, n,m < 0, 
{/?., /™} = 0,m<0,n>0orn<0,m>0. 

Let subspace M. S , P C g^ be defined as follows: 

p 

m=— s+1 

Brackets { , } could be correctly restricted to A4 S)P . It follows from the next 
Proposition: 

-p-l oo 

Proposition 3.2 Subspaces J p>s = Yl (dn)m + Y (flw)m are ^ea/s m g^. 

m=— oo rra=s 

Proof. It follows from the explicit form of commutation relations in g^. 
Now we are ready to prove the following important theorem: 



8 



Theorem 3.1 Functions {I^(L)} commutes with respect to the restriction of the 
brackets { , } on M. s>p . 

Proof. It follows from the combination of Kostant-Adler scheme and previous 
Proposition. Indeed, due to the fact, that {i^(L)} are Casimir functions on g^ 
they form a commutative subalgebra with respect to the brackets { , }o @]- Hence 
they will stay commutative after the restriction to M. S , P = (0«/i7p,s)*> due to the 
fact, that projection onto quotient algebra is a canonical homomorphism. 



Hamiltonian equations will have the following form: 

dl k 

-£ = {lb,H(®)}o, (20) 

where hamiltonian H is one of the functions or their linear combination. These 
equations could be written in the Lax form ||: 

^ = [L(w),M(w)), (21) 

where L(w) G A4 S , P , and second operator is defined as follows: M(w) = (P_ — 
P + )VH(L(w)). Here P± are projection operators on the subalgebra g^, 

vh(l(w))= £itm x « h (22) 

k=—p ij=l l 3 

is an algebra- valued gradient of H. 



4 Integrable systems in finite-dimensional quo- 
tients. 

The most interesting from the physical point of view examples usually arise in the 
spaces M. s ,p with small s and p. We will consider the case \s + p\ < 2. We will 
assume, that curve Tt is nondegenerated, i.e. a. t ^ a 3 - for i ^ j. This requirement is 
necessary for completeness of the family of constructed commuting functions. 

The basic algebra in all examples will be g = so(n), but analogous results are 
valid for g = gl(n) and g = sp(n). We chose for all examples so(n) algebra only 
because correspondent integrable systems are the most direct generalizations of clas- 
sical integrable systems connected with so(4) and e(3). 

4.0.1 Generalized interacting tops. 

Let us consider subspace jMy. In the case ai ^ 0, and as it follows from the 
arguments below, (g^/j7i,i) — g © g, hence A^i,i = (g © g)*. Corresponding Lax 
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operator L(w) G -M ,i has the following form: 



WiWj 



Let us again consider so(n) case and put = X^. Lie-Poisson brackets between 
the coordinate functions iff have the following form: 

X/(0) /(0) X _ r ,(0) , r ,(0) r ,(0) r ,(0) 

r/(!) ;(i)i _x /(!) a /W-lA /W A /W 

/ {0) 

Putting fej = a t 1//2 and making the change of variables: Uj = if], = we 

bibj 

obtain canonical coordinates of the direct sum of two algebras so(n): 

{mij,™>k,i} = Skjm^i - Sijirikj + S^m^i - S kii m jth 

{h,j, = dk,jk,i ~ + — $k,ilj,i, 
{hi, ma} = 0. 

Commuting integrals are constructed using expansions in the powers of w of the 
functions: Ik(w) = Tr(L(w)) h . As in the previous example, we are mainly interested 
in the quadratic integrals. Let 

n 

h( W ) EE I 2 (W) = hs(lff)w S = II ( W ' ^ l ff + Wlff) 2 . 

s=0 ij k^i,j 

By straightforward calculations we obtain: 

h = (-ir 2 t ^^(iffy 

i,j=l tt'lttj 

h x = jr ( y, aia2 "' an y^ > _ 2 a i Q 2--- a " ^(o)^(i) 

i,j=l k^i,j &k QjiQjj (lidj 



hn-i = " E (E a fc - (a, + aj )){lfff - *W 

i,j=l k=l 
hn=t (iff) 2 - 

Making the described above replacement of variables we will have: 
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n 

ho = (-ir 2 (&?&l- b 2 n ) E m% 

h = (-1)- 1 E ( E ^%^)K,) 2 - 2%^W. 



hn-x = ~ E (E b 2 k - (M + b)))l% - ■2l,,l, jWlJ l ll 

i,j=l k=l 

n 

It is evident, that functions ho and /i n are invariants. Functions hi,...,h n -i generate 
non-trivial hamiltonian flows. For the hamiltonian of the generalized interacting 
rigid bodies we can take either h n _i or hi. Correspondent M-operator and Lax 
equations are calculated straightforwardly. 



Steklov-Veselov (n = 3) case. In this case, making standard of variables ij 
m i — £ijk m ij we obtain the following set of commuting functions: 



^■ijk^ij i "' i ' ///.-" ' './ 

n n 

ho=Y, m l h 3 =J2 ll, 

k=l k=l 



E ( blb jf 3 )( m k) 2 - 2b k m k l k 



k=l 

n 



i,j=l Ofc 



Here ho, h 3 — invariant functions and ft,i,/i2 are two independent integrals discovered 
by Veselov [IB]. Commutation relations between coordinates l k and m k of so(4) are 
standard: 

{h, lj) = e ijkh, {m, raj} = e ijk m k , {m h l 3 ) = 0. 



4.0.2 Generalized Steklov-Ljapunov case. 

The last class of integrable systems we wish to consider here will be integrable 
systems that generalize classical integrable system of Steklov-Ljapunov on e(3) = 
so(3) + R 3 . Let us consider subspace A4o,2 = (Qn/^fi)*- It is easy to show that 
Mo,2 = (fl + fl)*- Corresponding Lax operator L{w) G M.0,2 has the following form: 

1,3=1 J 

We will again be concentrated on g = so{n) case and put X*- = Xij. Lie-Poisson 
brackets between coordinate functions are following: 

r/W AVX-X /« A /W_lA /W r /M 1 „ /(2) r ,(2) , r ,(2) r ,(2) 
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Commuting integrals are constructed using expansion in the powers of w of the func- 
tions: I k (w) = Tr(L(w)) k . We are again interested mainly in quadratic integrals. 
Let 

n 

h(w) = I 2 {w) = w 2 Y, W$V = w 2 II ( w - a ^f + w 2 - 

s=0 ij k^i,j 

Direct calculations give: 

h 2 = ± ^i^ a (D )2 



i,j=l di&j 



K + i = (-1)( E (E o* - (a* + ) 2 - 

i,j=l k=l 

n 

h n +2 = E ) • 

Change of variables: iff = Uj — 1/2 (a, + aj)pij, iff = p^ transform described above 
brackets to the canonical brackets on the half-direct sum so(n) + so(n): 

\hj, hi} = Skjhi — Suhj + djilki — Siklji, 

{kj,Pki} = hjPu - SiiPkj + SjiPki - SikPji, 
{Pij,Pki} = 0. 

After such transformation we obtain the following set of hamiltonians: 
h 2 = (-I)"" 2 E aia2 "' ^ {k - 1/2(0* + a^f 

i,j=l diOij 



h n +l = (-1)(E a k)( E Pij) ~ 2( E ^iiPv) 
n 

= E Pij- 

i,j=l 

Last two functions are invariant functions. First n — 1 give non-trivial flows on q*. 
We will chose function H = h 2 for the hamiltonian function. Correspondent M 
operator is: 

n 

M(w) = 2 22 — —{hj ~ l/2(ai + a j )p ij )w~ 1 w i w j Xi j . 

%,j=\ ai<x i 

Lax equation have the standard form ([[]). 
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Steklov-Ljapunov (n = 3) case. In this case, making standard replacement of 
variables Ij. = e^khj, Pk = tijkPij we obtain: 

n 

H = (-l) n ~ 2 a k(h - l/2(oi + a 2 + a 3 - a k )p k ) 2 . 

k=l 

Up to the rescaling of momenta pf Pi — > lapi hamiltonian H coincides with the 
hamiltonian of Steklov-Ljapunov system in the form of Kotter []15|. 
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